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I.  INTRODUCTION 


In  this  paper,  the  pressure  correlation  function  for  a single  boun- 
dary layer  over  a flat  plate  will  be  derived  for  the  subsonic  flow  regime. 
Numerical  results  will  be  given  based  on  experimental  measurements  of 
the  magnitude  of  the  mean  and  fluctuation  velocities  in  the  boundary  layer. 
One  application  for  this  work  is  that  of  the  propagation  of  a laser  beam 
through  a turbulent  boundary  layer.  Another  important  application  is 
that  of  the  generation  of  sound  by  boundary  layers  over  flat  surfaces. 

There  has  been  a tremendous  amount  of  work  directed  toward 
questions  relating  to  such  boundary  layers.  This  work  has  in  the  main 
been  experimental,  because  of  the  complexity  of  the  equations  governing 
fluid  flow.  Kraichnan  * ^ has  produced  theoretical  estimates  of  pres- 
sure correlations  at  the  boundary.  Other  literature  can  be  found  in 
standard  references. 

Index  of  refraction  fluctuations,  important  to  optical  propagation, 
axe  primarily  dependent  upon  density  changes.  Density  fluctuations  can  be 
caused  by  pressure  and/or  by  temperature  fluctuations.  As  a result  of 
recent  experimental  information  and  some  estimates,  it  has  seemed  rea- 
sonable to  concentrate  first  upon  the  effects  of  pressure  fluctuations.  For 
now  we  neglect  temperature  effects  and  treat  these  effects  in  a later  paper; 
in  effect  we  suppose  that  the  flow  is  adiabatic.  We  assume  air  is  an  ideal 
gas.  A sketch  of  the  boundary  layer  is  shown  in  Fig.  1. 
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Since  we  are  interested  here  in  subsonic  flow  regimes,  it  is  possible 
to  treat  the  fluid  as  approximately  incompressible.  It  is  true  that  the  mean 
flow  Mach  number  approaches  unity.  Nevertheless,  the  velocity  fluctuations 
(determining  the  pressure  fluctuations)  are  considerably  smaller  than  mean 
flow  values  and  their  Mach  number  is  small.  (The  question  of  whether  or  not 
one  can  treat  a flow  as  incompressible  is  governed  mainly  by  the  Mach  number.  ) 
The  incompressible  equations  of  motion,  the  Navier -Stokes  equations,  and  the 
equation  of  continuity  can  be  written 


and 


u + u • V u = 


1 2 
— p p + vV  u 


V . u = 0 


(1) 

(2) 


The  notation  in  these  equations  is  standard.  The  density  pQ  is  assumed  con- 
stant; the  dot  over  the  velocity  in  Eq.  (1)  represents  the  time  derivative  of 
that  quantity. 

The  equations  of  motion  are  of  course  very  complicated  at  this  level 
of  generality.  For  turbulent  boundary  layers  there  are  some  useful  and 
realistic  simplifying  assumptions.  The  boundary  layer  grows  very  slowly 
in  the  down-stream  direction:  the  characteristic  length  in  that  direction  is 
much  greater  than  the  thickness  of  the  boundary  layer.  To  good  approximation 
we  may  assume  that  the  averages  of  all  quantities  are  functions  only  of  the 
direction  perpendicular  to  the  wall;  the  mean  flow  in  particular  depends 
mainly  on  that  direction.  Consequently,  the  drag  on  the  flow  is  approximately 
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constant  in  the  down-stream  direction;  hence  the  pressure  gradient  (which 
overcomes  that  drag)  is  likewise  constant.  The  Poiseuille  flow  problem 
fulfills  these  approximations  exactly  and  is  discussed  in  this  connection  in 
Appendix  A.  We  assume  that  the  flow  is  statistically  stationary,  so  that 
averages  are  independent  of  the  time  at  which  they  are  taken. 

The  definitions  of  some  quantities  to  be  discussed  here  are  shown 
in  the  figure.  The  y-direction  is  perpendicular  to  the  plane  of  the  page. 

The  components  of  the  position  vector  and  of  the  velocity  in  the 
three  Cartesian  directions  are 

r_  = (x,  y,  z)  = (r  , r^,  r3) 

and  (3) 

u = ( u,  v,  w)  = (uj , u2,  u3)  ; 

the  choice  of  component  representation  is  made  as  convenient.  Other 
vectors  below  have  similar  representation.  We  write  the  velocity  and 
pressure  in  terms  of  mean  values  and  fluctuations 

u=u  + u/  = u ( z ) i + u7 

(4) 

p = p + p' 

Here^  is  taken  at  £ and  (later)  u’  at  r " and  similarly  for  p‘  and  p*. 

The  essential  method  to  be  employed  in  this  work  is  the  following: 
take  the  divergence  of  the  incompressible  Navier-Stokes  equation.  One 
finds  a Poisson  equation  for  the  pressure  where  the  right  tide  (source  term) 
under  our  assumptions  consists  of  second-order  velocity  fluctuations  and 
terms  ^uartic  in  such  fluctuations.  The  former  have  mean-velocity  gradients 
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as  coefficients.  The  Poisson  equation  is  solved  for  the  pressure  fluctuations 
the  correlations  (and  variances)  of  the  pressure  are  constructed  from  this 
solution.  The  quantities  so  found  are  integrated  numerically. 

In  Section  II  the  pressure  correlation  function  will  be  derived  with 
contributions  retained  through  the  fourth  order  in  velocity  fluctuations. 
Section  III  is  concerned  with  evaluation  of  the  fourth-order  contribution  to  the 
pressure  correlation  while  Section  IV  deals  with  the  second  order  term. 
Numerical  techniques  used  in  calculating  the  correlation  function  are  given  in 
Appendix  E.  Numerical  results  and  conclusions  are  presented  in 
Section  V. 


II.  DERIVATION  OF  THE  PRESSURE  CORRELATION 
From  the  above  assumptions  (see  Appendix  A) 


— • 2,  / . 

p = -ax-  w (z  ) 


(5) 


where  we  refer  pressures  to  the  pressure  at  the  origin.  Now  take  the 
divergence  of  Eq.  (1)  using  the  incompressibility  condition  giverr'by  Eq.  (2) 
and  find  (with  v'  2 the  Laplacian  in  r ' - space) 


v'  2p(r'  , t)  = -p 


u (r'  , t)  u (/  , t) 

o ' a P - 

Bra9r0 


(6) 


written  at  the  point  . We  use  the  summation  convention:  a repeated  index 
in  a single  term  is  to  be  summed  from  one  to  three.  If  now  we  substitute  Eq. 
(4)  in  Eq-  (6)  and  average  we  obtain  an  expression  for  the  average  pressure, 
noting  from  Eq.  (4)  that  averages  of  the  fluctuations,  n and  p ' , vanish. 
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Subtract  the  averaged  equation  from  Eg.  (6)  to  find  the  following  relation 
for  the  pressure  fluctuation  (the  time -dependence  is  implicit). 


^ P'  ' ‘ P°  ^ ‘-1  “ + % UB  - <<r  UV] 

<y  e 

HoB  <£> 


(7) 


3 - P 


o A a / 
dr  ct  8r  S 


Vfe  use  the  brackets,  < >,  or  the  overbar  for  time  averages,  as  convenient. 
The  last'  (average)  term  in  Eq.  (7)  is  a function  of  a'  alone  under  our  boundary 
layer  assumptions.  The  Kronecker  delta  has  values:  8y  . l,  i . j and  equal 
zero  otherwise. 


In  order  to  solve  Eq.  (7)  we  need  the  Green's  function  for 


the  problem. 


It  is  recalled  that  the  flow  is  confined  to  the  half-space  above  the  plane  (see 
Fig.  1).  The  boundary  conditions  on  the  pressure  at  the  plane  are  not 
Simple  and,  in  fact,  stated  exactly  would  involve  fluctuating  quantities.  An 
approach,  which  in  son.e  ways  simplifies  the  treatment,  is  the  following. 

First,  our  problem  consists  of  a plane  which  exert,  drag  upon  the  flow,  and 
the  compticated  flow  above  the  plane,  confined  of  course  mainly  to  the  boun- 
dary layer.  This  problem  can  be  replaced  by  an  equivalent  problem.  Suppose 
fha.  we  construct  a flow  below  the  plan,  which  is  the  mirror  image  of  the  actual 
one  above  it,  and  replace  our  problem  by  the  original  flow  plus  the  imaged 
flow  plus  a force;  the  boundary-layer  plane  of  course  exert,  a drag  upon  the 


flow  and  when  w 
Fortunate  ly  this 


e remove  that  plane  we  must  insert  a force  equal  to  that  drag, 
force  is  negligible  (see  Appendix  B).  We  proceed  with  the 
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situation  as  described  consisting  of  the  actual,  physical  flow  and  the  mirror- 

imaged  flow.  Then  the  Green's  function  needed  for  the  solution  of  Eq.  (?)  is 

the  free  space  Green's  function  for  the  Laplacian  operator--essentially 

Further,  we  suppose  that  all  fluctuating  quantities,  pressure  and  velocity, 

vanish  sufficiently  far  outside  the  boundary  layer  regions,  whether  in  the 

physical  flow  or  the  mirrored  flow.  Then  all  solutions  of  the  homogeneous 

equation  derived  from  Eq.  (7)  (that  is,  solutions  of  Eq.  (7)  with  the  right 

4 

side  equal  to  zero)  must  vanish.  The  solution  of  that  equation  becomes 


H. 


/ O'  0 

ar0tt  3r0B 


(r/  ) dr' 
"0  “0 


with , as  above, 

h,b  = 2“(z'> 6„  1 u'b  + u'<*  u's  - <u'« “V 


(8) 


The  velocity  functions,  and  integrals  involving  them,  extend  now  over  all 
space  (they  are,  as  stated,  primarily  confined  within  the  boundary  layer 
and  the  mirrored  boundary  layer). 

To  find  the  pressure  correlation  from  this,  write  Eq.  (8)  for  the 
pressure  at  the  point  r_"  and  multiply  that  equation  by  Eq.  (8)  at  the  point  r ' 
and  average  to  find 


Cpp=  <p'(r')P  (£  >> 


L - Lo' 


-1 


(9) 


5r0u  Ss  SrOvir06 


Hjt'o1  Hv«  d£od£o 


The  question  of  convergence  of  these  integrals  is  of  primary  importance. 

H 

Suppose  the  integral  over  is  performed  first;  we  can  expect  that  the  cor- 
relation in  the  integrand  will  vanish  for  sufficiently  different  from  r^, 
so  that  the  integral  will  be  convergent  at  infinity.  Furthermore  there  is  a 
question  of  convergence  when  r_  - r ^ vanishes.  The  differential 
element  of  volume  d£g  will  take  care  of  this  apparent  singularity.  There  is 
however,  more  difficulty  with  the  outer  integral  over  r^ , in  particular  at 
infinity  for  that  integral.  This  is  a familiar  problem  if  one  attempts  to 
deal  seriously  with  a pressure  representation  as  given  in  Eq.  (8).  The  fact  is 
that  there  is  a considerable  amount  of  cancellation  in  that  integration  process 
which  results  from  the  fact  that  the  integrand  of  Eq.  (8)  is  a quadrupole  source 
for  the  pressure  field.  The  proper  way  to  deal  with  this  difficulty  is  to  inte- 
grate Eq.(9)  by  parts  using  the  derivatives  with  respect  tc  r^  . This  is 
what  will  be  done  below.  When  that  manipulation  is  carried  out  we  have 
the  factor  | t/  - r^j  " replaced, for  r^  large, by  something  like  |r^|~ 
which  will  give  a convergent  result  at  infinity.  Then  one  might  expect 
difficulties  at  the  singular  point  |jr^  - r^j  = 0.  Fortunately,  the  integrcJ 
is  also  convergent  there,  because  of  angular  cancellations.  The  difficulties, 
which  arise  at  the  singular  point  when  we  integrate  by  parts,  are  discussed 
in  Appendix  D. 
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Consider  the  correlation  in  the  integral  (9).  It  consists  of  second- 
order,  third-order  and  fourth-order  velocity  correlations,  the  second-  and 
third-order  with  coefficients  involving  the  mean  flow.  To  begin,  we  note 
that  third-order  velocity  correlations  are  usually  very  small  compared  with 
even  order  correlations.  Indeed,  for  a velocity  field  with  Gaussian  statistics, 
odd  moments  would  vanish.  In  actual  flows  the  odd  moments  of  fluctuating 
quantities  are  typically  very  small.  We  neglect  them. 

There  have  been  many  measurements  made  of  second-order  velocity 
correlations  but  rather  few  of  fourth-order  correlations.  We  make  the 
reasonable  assumption  that  the  velocity  field  is  not  too  far  from  Gaussian; 
then  the  fourth -order  velocity  correlation  can  be  written  as  follows 


<uou V Vfl > = Q*V (^0'  4>  qS6 

+ Q(v6(~o'  -o) 

+ Qa0  (-O'  -0^Qy(>  {-O'  ^0)  ’ 


(10) 


the  result  of  this  Gaussianity  assumption  is  that  the  fourth-order  velocity 
correlation  can  be  written  in  terms  of  quadratic  expressions  involving  the 
second-order  correlation.  We  have  defined  the  second- order  velocity 
correlation: 


We  let  the  pressure  correlation  be  represented  by 


PP 


= c*1*  +c{2* 


PP 


PP 


(11) 


(12) 


where 


r(l)  /Po  V ff,  i /.-l  * * 

pp  =l^7  Jjl-  ' -o ! ±'  r-o 


dro»*roe  aroYdro6 


X [4J|i'0)u|i0)Jffl  SlSfl^O'io)]  d£o  d4 

r (2>  _ ( p°\z  ff,  ' / . “ 1 1 # * -i  a4 

PP  U"/  _£ol  ^ _£ol  &r'  Sr'a&r'  hr* 


(13) 


rO»  r0pwr0v  06 


(14) 


x [2Ryv<i'o-  4>  Q96(i'o  • i>] d4  d4  • 


efficients.  One  of  his  conclusions  was  that  the  first  integral,  C ^ 

PP 


In  this  last  equation  the  symmetry  under  index  interchange  has  been 

utilized.  Further,  the  last  term  in  the  fourth-order  correlation  given  in 

Eq.  (10)  (that  last  term  representing  simple  averages)  is  cancelled  by  the 

1-2 

last  term  in  the  function  given  in  Eq.  (8).  Kraichnan  has  con- 

sidered the  problem  of  the  pressure  correlation  at  the  boundary -layer  wall. 

He  made  some  progress  in  the  estimation  of  the  size  of  the  correlation  co- 

\ was 

more  important  than  the  second  discussed  here,  the  reason  being  that  the 
first  involves  the  mean-flow  velocity  whereas  the  second  involves  only  the 
(smaller)  velocity  fluctuations.  We  carry  both  of  the  contributions  to  the 
pressure  correlation  and  will  determine  through  computation  whether  the 
assertion  is  justified  or  not.  There  is  some  reason  to  question  it,  despite 
its  apparent  reasonableness.  Measurements  of  rms  pressure  fluctuations 
at  the  wall  give  values  which  appear  to  be  j pQ  u#  rather  than  ^ pQ  u*u , 
as  would  be  the  case  if  the  mean  flow  term  dominated  the  result.  To  proceed, 
further  (reasonable)  simplifying  assumptions  are  required.  First,  for  the 
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L 


integral  c^^  we  suppose  that  the  second-order  velocity  correlation  can  be 
written. 


/ » 


with 


°B6  (-0’-0)  = u (z0)  u (zo)  °B6 

io  = i = 4 • 


«W0>  = ( T H.#2  (Zft)) 


1/2 


and 


",  " V /I  »Z  . • . V 

u(z0)  = <ju  (zQ)  ) 


1/2 


(15) 


jje  t 

where  Q^gtr^)  is  the  normalized  second-order  velocity  correlation  for  homo- 
geneous and  isotropic  turbulence  (normalized  by  the  rms  values  of  the  velocity 
fluctuation  in  one  direction).  Here  we  have  used  a normalization  involving  the 
velocity  fluctuations  at  the  two  points  and  rj.  The  normalized  second-order 
velocity  correlation  has  an  outer  scale  which  is  approximately  the  distance 
from  the  wall;  for  that  scale  we  choose  M.  The  scale  M is  defined  as 
M 1 = [(z'0)  1 + T 1]  1 where  T is  the  thickness  of  the  fluctuating  boundary 

layer,  defined  as  the  distance  from  the  wall  where  u attains  0.  9 of  its 
maximum  value.  This  T is  greater  than  6 (the  displacement  thickness  defined 
ky  ® ^ ^ * u ^ was  felt  that  this  value  T was  appropriate  here, 

rather  than  the  smaller  6,  The  substance  of  this  assumed  form  is  that  we 


propose  that  the  second-order  velocity  correlation  is  locally  isotropic,  but 
corrected  for  the  variation  in  strength  of  the  velocity  fluctuation  over  the  dis- 
placement vector  r_Q.  The  nature  of  the  second-order  velocity  correlation  for 
homogeneous  and  isotropic  turbulence  is  well  known  through  wind  tunnel 
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measurements.  Later  we  shall  take  numerical  values  from  standard 
literature  references  ^ . In  fact,  for  our  incompressible  flow  problem, 
the  correlation  function  becomes,  for  large  Reynold's  numbers,  essen- 
tially independent  of  the  Reynold's  number:  a universal  function.  Con- 
sequently, for  different  boundary  layer  measurements  it  will  not  be 
necessary  to  change  these  functions.  The  normalized  velocity  correlations 

c 

for  homogeneous  and  isotropic  turbulence  can  be  written 


Q. . (r ) = 

ij 


= f(r/M)  - g(r/M) 


r.r . + g(r/M)6. . 
i J * ij 


(16) 


The  f and  g and  Q. . (r)  are  simply  behaved  functions  of  r/M  where  M,  the 
^ J 

"distance"  from  the  wall  is  the  outer  scale  for  the  velocity  correlation  as 
proposed  above. 

For  the  second-order  velocity  correlation  (dropping  the  scaling 
factor  for  now)  it  will  be  convenient  sometimes  to  use  the  notation 

Qy  <£.)  = Qj(r)r.rj  + Q2(r)6„  (17) 


with  the  functions  Q^,  f and  g related  as  shown.  Finally  for 
incompressible  flow  one  can  write  the  function  g in  terms  of  the  function 

5 

f as  follows 


g(r)  = f + 1/2  r f' 


(18) 


and  of  course  there  is  a corresponding  relation  between  and  Q2.  The 


relation  Eq.  (18)  follows,  as  indicated,  directly  from  incompressibility. 
That  condition  can  be  written 


3 « * 3 ❖ 

57“  Qij  * ST-  Qij  = 0: 


(19) 


The  relation,  Eq.  (18)  is  obtained  from  this  by  differentiating  Eq.  (16). 
If  we  substitute  from  Eq.  (15)  to  Eq.  (19)  into  Eq.  (13)  we  find 


r(l)  /pQ ff  , ' , » • / -I  * » 

cpp  -(15T7  )Jd- od£o  |i  -io|  |i  -£o 


-1 


~ U (Zo)Ju'(z^ 


L dz 


3 — . » . 

u (zo) 


0 


dx 


7T  Q33  0 “ — 0 


",  * s 

u <v 


(20) 


In  order  to  proceed  further  with  this  term,  simplifying  assumptions 
must  be  made.  For  now  we  will  consider  the  more  complicated  fourth -order 
contribution  and  return  to  the  evaluation  of  Eq.  (20)  in  Section  IV. 

III.  DETERMINATION  OF  THE  FOURTH -ORDER  CONTRIBUTION 
TO  THE  PRESSURE  CORRELATION 

(2) 

Continue  now  with  a discussion  of  the  term  C given  in  Eq,  (14).  For 

PP 

(2) 

the  term  C we  take  for  the  second-order  velocity  correlation  the  following 


2 

where  our  boundary  layer  assumptions  mean  that  u c is  a function  of  only 
z*  . Substitute  this  relation  for  the  velocity  correlation  into  Eq.  (14).  Then 
integrate  by  parts  and  use  Gauss'  theorem  to  obtain  the  result  (see  below  for 


I), 

<22> 

From  Gauss'  theorem  one  also  obtains  integrals  over  the  surfaces  bounding 
the  region  of  turbulence.  These  surfaces  are  taken  outside  the  physical  flow 
and  the  image  flow;  on  those  surfaces  all  the  fluctuations  vanish  so  the  sur- 
face integrals  vanish.  There  is  also  the  question  of  the  singular  points  which 
arise  when  the  denominators  of  Eq.  (10)  vanish.  It  is  shown  in  Appendix  D 
that  these  singularities  in  fact  cause  no  difficulty.  We  have  used  new  variables 
in  Eq.  (22)  (as  above  and  as  follows), 


//  / 

Lq  ' Lq 


and 


1=  £.  - £0 

and  defined  the  integral 


(23) 


' fA to  I 1 - lo  I 'lQav<^)  QeVlo>- 


(24) 


Substituting  from  Eq.  (18) into  Eq.  124)  we  identify  the  following  different  types 
of  terms  for  that  integral: 


I , I<a>  +I(bl>  +I(b2)  +I(C) 

oyB6  ovBS  tyySfi  <YyB6  oySS  * 


r0o  r0B  r0y  r06  ' 


‘/oXoli-Iol  ’*  Ql2  <r0> 

= A1<5)5»!6!v;6+A2(5)  [5aVv»  + ] 

+ A3(!)  ^6a6  «yS  + 6b#  + iai  6ay  j, 

where  the  coefficient  of  A£(§)  has  the  six  permutations  of  the  indices,  two- 
by-two 

4yB6  = J d-^0  U " -0  I Ql(r0)  Q2  (r0)  r0a  r0y  606 

= [Bi<5)^  VB2(?>6ov]6B5  ; 


(25) 

(26) 

(26a) 


r (b2  ) = (bl)  . 
ovBS  BSo-v 

XjyB6  J d^-0  | ~ ~ ^-0  | Q2  (r0}  6av  6B6 


= C (")  6 6b.  . 

QfV  Bo 


(27) 
(27a) 

(28) 
(29) 
(29a) 


We  know  the  integrals  must  be  of  the  forms  of  Eqs.  (26a),  (27a)  and  (29a) 
because  of  their  vector  (transformation)  properties.  One  must  now  find  the 
functions  Ay  A^,  Ay  By  and  C.  These  functions  are  evaluated  in 
Appendix  C.  For  substitution  in  the  integrand  of  Eq.  (22)  we  must  evaluate 
the  second-order  derivatives  of  the  integrals  given  in  Eqs.  (25)  - (29a). 
Again,  by  the  vector  properties  of  the  expression  we  know 


wiwwn 


m i 


■Trv"S?;"  IayB6(-)  = D1<?>?*S  + D2<5>6«|3  ' 

Substituting  the  results  given  in  Appendix  C into  Eq.  (30)  and  adding  the 
coefficients  of  corresponding  terms,  we  find  for  the  functions  and 

Dj(§)  = 4tt  Ex(?)  + | ?'5  E6  + 2§'5  F4  + 3§“5  G2 


(30) 


- j’?Q1(§)+  ?_1Q2(?)]  [ 

D2(5)  = 4nUs2El  + | E3  - £ I'3  E6  + | F, 


(31) 


— P”^F  - C G 

3 * *4  5 


These  expressions  may  be  simplified  by  the  use  of  Eqs.  (16)  - (18) 
yie  Iding 

Di<?>=  ^ |12  Ei<5>+  7?“5e6(§)| 

2 

d2(?)=  |l2E1(?)+  7?  5 E6(§)+  10r2E3(5)+  10?~2f2(§)J 


(32) 


Ola) 


(32a) 


Equation  (31a)  unlike  Eq.  (31)  is  easy  to  evaluate  at  5 = 0 and  differentiation 

w.  r.  t.  5 is  more  easily  accomplished. 

Here  E (x),  F (x)  and  G (x)  are  defined  by  the  expressions 
n n n 
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rn  q{r)  dr,  n odd 


and 


(33) 


A 

J 


r11  q(r)  dr,  n even 


with  q(r)  = Q2(r)  for  En;  q(r)  = Q^r)  Q2(r)  for  Fr;  q(r)  = Q^r)  for  Gn. 
Substitute  Eq.  (30)  into  Eq.  (22),  use  the  known  relation 


2 |r'  -4*6  (r'  - j^), 


(34) 


define 


c(2)  _ c(2a)  + c(2b)  + C(2c) 
PP  PP  PP  PP 


(35) 


and  find  for  the  term  involving  D2> 

, 2 


\ 4 

J 2 [u'  (z'  ) ] D. 


(r) 


(35a) 


with  r_=  r"  ~ L • The  other  contribution,  from  Dj,  JCpp^  and  CppC^J  is 


■ (2c)"l 


more  complicated.  To  simplify,  use  the  definition  of  5 [Eq.  (23)J  to  write 


-*2  -1  2 
r^  - £ = r_  - 5.  and  noting  that  — - — — | £ - | = ^ | r_  - §_|  find 

dr_9„  a *5 


-1  3[£  . (r  - S)]2  ^ 


r - II 


r -II 


(36) 


Using  5 • r = j (?2  + r2  - |^r  - Sj  2)  we  have  for  the  right  hand  side  of  Eq.  (36), 
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(37) 


3 

4 


/ 2 rV 

(I  -s  ) 


Li£  " ? 


Substitute  Eq.  (37)  and  the  first  term  of  Eq.  (30)  into  Eq.  (22)  to  find  (aside 

from  the  delta  function  contribution  at  the  singular  point  r*  - rl,  given  below) 

— 0 


where  £_=_£+  Jk  and  M - (T  ~ * + | z " - gj  j ~^)  "^.  Note  the  scaling  factor  used 
in  Eq.  ( 1 6)  has  been  used  explicitly  in  Eq.  (38).  This  scaling  also  applies  to 
Eq.  (35a)  and  Eq.  (38b)  though  it  is  not  shown  there.  The  apparent  singularity 
in  the  neighborhood  of  _r  = in  the  jj  and  5^  integration  is  in  fact  in  integrable. 
This  is  most  easily  shown  by  using  spherical  coordinates  about  the  origin; 
convergence  results  from  the  angular  integrations. 

The  last  contribution  comes  from  the  delta  -function  at  _r  = 5;  care  must 
be  exercised  in  order  not  to  lose  this  portion- -as  cautioned  at  the  end  of 
Appendix  D.  It  can  be  found  by  taking  the  integral  of  the  term  involving  D^(g) 
in  the  volume  |r/  - r^  | s e,  e very  small.  We  have  for  it  from  Eq.  (22) 

f and  the  first  term  of  Eq.  (30), 

I: 

y. 


c‘2c>  = f 

pp  \41r/  J. 

\ / I / / 


dio  u/4(zo) 


r -^1  * e 


a 


ar'  ar' 

a 0 


li  -£i 


Dl<!>  tah-  (38b) 


Following  the  result  given  in  (D-7)  this  is 


„(2c)  po2  / 4.  / 2 , . 

C'  = - -rr  u z ) r D.  (r). 
pp  Off  I 


(38c) 


In  Appendix  F,  the  techniques  used  in  the  calculation  of  Eq. 

(38)  will  be  discussed.  Consider  now  the  second  order  contribution. 

IV.  DETERMINATION  OF  THE  SECOND-ORDER 
VELOCITY  CONTRIBUTION  TO  THE  PRESSURE  CORRELATION 

It  is  plausible  to  expect  that  the  second-order  velocity  fluctuation 
terms  will  dominate  the  fourth-order;  this  has  often  been  suggested  in 
the  literature  (see  e.  g.  Ref.  1).  We  shall  see  from  our  computations  that 
this  is  usually- -though  not  always- -the  case. 

In  the  determination  of  the  second-order  contribution  several  different 
assumptions  were  considered.  The  assumption  finally  chosen  as  the  most 
tractable  was  that  the  prime  location  replaces  the  double  prime  location  except  ir 
Q 


33- 


- (*o> 


a*, 


» , a 

U (z0)-£. 

dx 


-l')  * 


*Z  33  0 --0 

0 


■■■■» ^ ii 


az 


/ u(*0) 


u'(z') 


a x 


72  Q33  ^ ^ 


(39) 
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This  is  reasonable  for  Q,_-*0  as  r.  - ri  becomes  large;  it  is  consistant 

53  — U — U 

with  the  approximation  made  for  the  fourth  order  contribution  in  Eq.  (21). 
Inserting  Eqs.  (39),  (23)  into  Eq.  (20)  yields 


d-Q  Q 

I ” 


33^ 


±0 


(40) 


Analogous  to  treatment  of  the  fourth  order  contribution,  define 


W5’ 


/ <1X0  1 1 - £o  I Q*e  <£o>  = Al'5>  V B + Aa<5»6a  B+  B2(5,6a  B 


(41) 


where  is  determined  first  by  using  the  double  cross-product  with  the 
vector  j[  and  then  Aj  is  determined  by  taking  the  inner  product  (summing 
on  the  indices).  Using  results  from  Appendix  C. 

B2(§)  = 4jt  | s'1  F2  + Fj  | 

A2<:>=f?  l5te'‘E4  + e3]  - + Ej | 

Al<5>  =1*  lt‘5E6  *^,1 

where  E,  E are  defined  as  in  Eq.  (33)  except  that  here  q(r)  = Q^(r)  for  E?n  and 


-22- 


q(r)  = Q,(r)  for  F . Using  these  expressions  and  Eqs.  ( 1 6)  - (18)  we 
c,  n 

have  the  scalar  function 


Substituting  Eq.  (42)  into  Eq.  (40)  and  taking  into  account  scaling  results  in 


cpp=”2  / —«,h  .'d  II*  / d * 

(43) 


where  0 is  the  angle  between  the  projections  of  and  r_  in  the  x-y  plane. 

2 2 2 

Noting  that  ^ = p cos  (a-  if)  with  ais  the  angle  between  the  projection 
of  _r  in  the  x-y  plane  and  the  x axis  unit  vector  allows  Eq.  (43)  to  be  written 
as 


C(1)(x, 

PP 


z') 


C(1)  (-), 
PP  1 


0,  i,  z') 


2 

cos  a 


C<>> <0, 
PP 


r\> 


i#) 


. 2 
sin  a 


(44) 


/ 2 2 

where  ri  = *x  + y and 


C C Lt 

is  replaced  by  p cos  $ for  the  argument  y 
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and  p sin  <j,  for  x = 0.  The  cross  term  vanishes  for  integration  over  $. 

V.  NUMERICAL  RESULTS  AND  CONCLUSIONS 
This  section  contains  numerical  results  and  conclusions  based 
on  calculating  the  pressure  variance  and  correlation  function  from  Eq. 

(38)  and  Eq.  (44).  For  the  interested  reader  the  details  of  the  calcula- 
tion may  be  found  in  Appendix  E and  supplemental  numerical 
results  may  be  found  in  Appendix  F.  In  order  tc  present  re- 
sults which  are  cohesive  a brief  step-by-step  procedure  for  obtaining 
the  numerical  results  is  presented  with  equation  and  figure  numbers 
prefixed  by  letters  referring  to  the  corresponding  appendix.  Some  of  the 
results  presented  here  should  more  reasonably  have  been  presented  in 
Appendix  F;  however,  they  are  presented  here  to  Bhow  the  reader  how  to 
calculate  the  pressure  variance  for  his  own  boundary  layer  data  without 
referring  to  the  cited  appendices. 

In  order  to  obtain  numerical  results,  all  intermediate  functions 
were  calculated  and  stored  in  the  computer  in  tabular  form.  Numerical 
integrations  were  preformed  by  Gaussian  Quadrature  integration  techniques 

with  the  interval  step  size  determined  by  an  accuracy  criteria  of  one  part 

3 

in  10  . Values  of  the  functions  not  in  the  table  were  obtained  by  linear 
interpolation  of  the  log  of  the  function  for  positive  functions.  Close  to  the 
origin  or  for  negative  functions  second  to  fourth-order  Lagrange  inter- 
polation techniques  were  used.  The  calculation  starts  with  the  definition 
of  the  longitudinal,  f(r),  and  transverse,  g(r),  velocity  correlations  (16). 
These  may  be  taken  from  experiments.  In  our  calculations  we  must 
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we 


differentiate  these  functions:  a smooth  representation  is  needed.  Thus 
use  the  calculation  which  has  been  made  using  the  Weiner  Hermite 
expansion^*.  Those  results  fitted  large  Reynolds'  number  experiments 
very  well  and  were  differentiable.  The  second  step  in  the  calculation 
is  the  definition  of  Qj(r)  Eqs.  ( 1 6) -( 1 7).  Using  f,  g,  and  the  functions 
Dj(r)  Eq.  (31a)  and  D£(r)  Eq.  (32a)  are  determined  by  numerical  integration 
Figure  2 shows  the  results  for  D^(r).  After  the  angular  integration  in  Eq. 
(38a)  has  been  performed  analytically  and  the  results  expressed  in  terms 
of  the  complete  elliptic  integrals,  a change  of  variables  and  an  integration 
by  parts  is  made  on  the  p coordinate.  This  technique  removes  a numerical 
convergence  problem  but  requires  the  evaluation  of  the  first  derivation  of 
Dj(r)  which  can  be  obtained  from  Eq.  (3la).  Once  the  flow  field  turbulent 
velocity,  u/(z),  is  determined,  the  fourth-order  contribution  to  the 
pressure  correlation  and  variance  can  be  obtained  directly  (see 
Appendix  E)*  Results  for  this  calculation  are  presented  later  in  this 
section.  As  an  aside  the  fourth-order  contribution  to  the  pressure  vari- 
ance can  be  written  as 

2 

C<2)  (z\  o)  = ~~  -[u'(z')]  O (o)  + 

pp  c.  Jr  o 

where 

M(xx)  = 


dC3  IV  (|z/  + g.jl  ] M(xx)M"  (E3) 
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f 

l 

k 

j XX  = ?3  M'1  (E4) 

i 

• I 

M = (T-1  + | z'  + §3  I"1)'1 

«► 

D2(0)  = -0,  54154 

and  T is  the  thickness  of  the  boundary  layer.  The  function  M is  plotted 
in  Fig.  3;  the  reader  may  obtain  the  fourth-order  variance  for  other, 
particular,  experimental  values  of  u,(z). 

We  now  consider  the  second-order  contribution  to  the  pressure  vari- 
ance and  correlation  functions.  From  Eqs.  (42)-(44)  it  is  seen  that  D(§) 

V 

is  required.  This  scalar  function  is  obtained  using  Eq.  (42)  and  the  defini- 

{ * ^ *7 

tion  of  S E^(?)  [defined  after  Eq.  (41)].  The  interested  reader  may  then 

\ obtain  the  second-order  pressure  variance  directly  from  Eq.  (43).  First  it 

f 

is  noted  that  £ = 0 implies  that  the  angular  integration  can  be  performed 

2 

analytically  yielding  a factor  n for  terms  not  including  5,  and  n/2  for 

1 A 

| 2 

terms  with  ^ . The  integral  over  p is  then  calculated  as  a function  of  ? 

! 2-1 

and  normalized  by  4tt  m . This  function  is  tabulated  in  the  computer  and 

used  in  Eq.  (43)  to  obtain  the  variance  for  any  boundary  layer.  This  func- 
tion is  plotted  in  Fig.  4,  In  order  to  obtain  the  second-order  correlation 
the  angular  integrations  for  Eqs.  (43)-(44)  are  also  performed  analytically 

and  lead  to  the  tabular  function  B(x)  Eq.  (E10).  The  calculation  of 
■ rr 

[ . then  proceeds  via  a double  numerical,  integration  (E12),  given  the  flow 
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field  data  u'  u'.  Graphical  results  for  f(r),  g(r),  Q^r),  D2(r),  Dj(r)* 
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and  § E^(§)  are  *ount3  *n  Appendix  F. 

In  order  to  obtain  specific  numerical  results,  flow  field  quanti- 
ties must  be  given. 

The  flow  field  quantities  necessary  to  perform  the  calculations 
are  the  normalized  mean  turbulent  velocity  used  for  the  calculation  of  the 
cjbp^  term  and  the  burbulent  fluctuation  velocity  in  the  z direction  multi- 
plied by  the  derivative  of  the  flow  velocity  in  the  x-direction,  which  is 

used  in  the  calculation  of  the  C ^ term. 

PP 

Our  calculations  were  carried  out  for  two  separate  boundary 

layer  flows:  first,  for  a thickened  turbulent  boundary  layer  relating  to 

7 

experiments  carried  out  at  NASA -Ames  . The  thickened 

layer  was  generated  by  a flow  of  approximately  Mach  0.  9 and  a Reynolds 
number  of  3 x 10^/ft.  The  thickening  came  as  a result  of  upstream  pins 
and  spoilers  projecting  from  the  wall  over  which  the  boundary  layer 
formed.  The  second  boundary  layer  for  which  calculations  were  made 
was  an  idealized  one  (called  here  the  canonical  boundary  layer).  Data 
were  taken  from  earlier  experimental  work  with  such  boundary  layers 
(Hinze,  Ref.  8>  Figs,  7-4  and  7-10).  The  flow  field  quantities  neces- 
sary for  the  calculation  are  presented  in  Figs.  5-8  vs.  the  normalized 
distance  from  the  wall.  The  normalizing  factor  for  all  distance  measure- 
ments is  the  displacement  thickness  6 defined  previously.  Figures  5-6  are 
for  the  experimental  case  while  Figs.  7-8  are  for  the  canonical  case.  For 
these  two  flow  fields  the  normalized  (normalization  factor  pQ  u /4) 
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Fig.  8.  U3  Turbulent  Velocity  * UB'  - Canonical  Boundary  Layer 


variance  and  correlation  functions  are  calculated.  In  Figure  9»  the 
normalized  pressure  variance  for  the  thickened  (experimental)  boun- 
dary layer  is  shown.  The  variance  at  the  wall  and  at  the  peak  (which 
here  occurs  at  the  distance  of  about  2.  25  displacement  boundary  layer 
thicknesnes)  are  larger  than  the  corresponding  values  (see  below) 
relating  to  the  canonical  boundary  layer,  probably  because  of  the  greater 
value  of  the  turbulent  fluctuations  farther  from  the  wall.  The  position 
of  the  peak  of  the  variance  is  considerably  farther  from  the  wall  for 
this  thickened  experimental  boundary  layer  than  is  the  case  for  the 
canonical  layer.  This  is  probably  again  a result  of  the  fact  that  the 
thickened  layer  shows  greater  turbulent  fluctuations  farther  from  the 
wall.  It  is  noted  in  Fig.  9 that  the  fourth -order  velocity -fluctuation 
contribution  to  the  pressure  variance  is  slight,  amounting  to  at  most  but 
a fraction  when  compared  with  the  second-order  contribution.  In  Fig.  10 
we  show  the  corresponding  variance  for  the  canonical  layer.  The  variance 
is  essentially  the  same  as  the  measured  values  given  in  Ref.  7.  Here  the 
fourth-order  contribution  is  a considerably  larger  portion  of  the  total 
pressure  variance  than  was  the  case  with  the  thickened  boundary  layer. 

The  suggestion  is  that  for  well  developed  idealized  boundary  layers  it  may 
not  be  valid  to  assume  that  the  fourth-order  contributions  may  be  neglected. 
It  is  further  noted  that  in  this,  canonical,  case  shown  in  Fig.  10,  the  peak 
of  the  pressure  variance  occurs  extremely  close  to  the  wall  (0.04  6), 
essentially  in  the  vicinity  of  the  laminar  sub-layer.  This  comes  about 
because  the  velocity  - fluctuations  tend  to  peak  sharply  within  that  layer.** 
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The  normalized  correlation  functions  are  shown  in  Figs. 

11-13  for  the  experimental,  thickened  boundary  layer.  In  Fig.  11,  the 
correlation  in  the  z-direction  normal  to  the  wall  is  shown  for  a base 
position  located  at  the  wall  and  in  Fig.  12  the  corresponding  correlation 
is  shown  for  a base  position  located  at  the  position  of  the  peak  variance 
for  that  boundary  layer.  In  both  cases  the  fourth-order  velocity  contribu- 
tion to  the  pressure  correlation  is  slight.  It  is  noted  that  the  correlations 
nave  a long  tail,  retaining  10%  of  their  value,  even  out  to  9 boundary  layer 
thicknesses.  In  Fig.  13  the  x-  and  y-correlation  functions,  thftt  is  in  the 
streamwise  and  transverse  (or  spanwise)  directions  are  shown  for  the 
experimental,  thickened  boundary  layer.  It  is  interesting  to  note  that 
the  x-correlation  shows  a marked  negative  lobe;  the  other  correlations 
remain  positive.  The  scale  for  the  three  correlations  is  essentially  the 
same.  In  Figs,  14-16  the  normalized  correlation  functions  are  shown 
for  the  canonical  layer.  It  is  seen  in  those  figures  that  (again)  the  fourth- 
order  velocity -fluctuation  contribution  to  the  pressure  correlation,  though 
smaller  than  the  second  order,  is  in  some  ranges  appreciable.  Scales  of 
the  correlation  functions  in  the  three  directions  are  essentially  the  same, 
the  displacement  thickness.  The  x-correlation  (in  the  streamwise  direc- 
tion) shows  a slight  tendency  to  go  negative,  but  this  is  not  viewed  as  relia- 
ble. It  is  interesting  that  the  measured  pressure  correlations  at  the  wall 
likewise  show  no  negative  lobe  (as  is  found  experimentally,  see  Willmarth 
and  Woolridge,  Ref.  9). 

To  sum  up,  the  fourth-order  contribution  to  the  variances  and 
correlations  appears  to  be  negligible  for  the  thickened  experimental 
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Fig.  12.  Normalised  Pressure  Correlation  in  the  z Direction  at 
2.25  Del  from  the  Wall  - Experimental  Boundary  Layer 
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turbulent  boundary  layer,  though  in  certain  regions  it  would  be  necessary 
to  include  it  for  the  canonical  boundary  layer  (where  the  velocity  fluctua- 
tions are  confined  for  the  most  part  to  regions  very  close  to  the  wall). 
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Appendix  A 


POISEUILLE  FLOW  PROBLEM 


There  is  a problem  which  exactly  satisfies  the  assumptions  made 
for  the  single  boundary  layer  problem.  It  is  the  flow  between  parallel 
plates  --  called  Poiseuille  flow.  Consider  the  turbulent  flow  between  such 
plates,  where  the  Reynolds'  number  is  large  enough  to  generate  two  (tur- 
bulent) boundary  layers,  one  associated  with  each  of  the  plates,  and  a 
central  region  of  lesser  fluctuation  intensity  between  the  plates.  It  is 
plausible  to  suppose  that  each  one  of  the  boundary  layers  is  related  to  a 
single  boundary  layer  over  a single  flat  plate.  There  are  some  differences 
between  the  two  fluid  flow  problems;  e.  g. , in  the  flow  between  flat  plates 
the  central  region  shows  a residual  level  of  turbulent  activity  which  exceeds 
that  observed  outside  a single  boundary  layer.  Further,  we  know  that 
a single  layer  increases  (slowly)  in  thickness  downstream;  this  doesn't 
occur  for  the  Poiseuille  problem,  in  equilibrium. 

In  our  application  we  can,  as  stated,  treat  the  fluid  flow  as  one 
which  is  statistically  stationary,  that  is,  a flow  in  which  all  of  the  (time) 
averages  are  independent  of  the  time  at  which  the  average  is  taken. 
Furthermore,  for  the  parallel  plate  problem  after  the  flow  has  developed, 
the  flow  may  be  assumed  to  be  statistically  homogeneous  in  planes  parallel 
to  one  of  the  plates.  All  averages  depend  only  on  z,  see  Fig.  A-l.  (It  is 
not  assumed  that  the  flow  is  statistically  isotropic.)  It  will  be  seen  that 
these  conditions  permit  considerable  simplification  in  the  statistical 
quantities  which  we  need.  The  geometry  of  the  flow  problem  is  sknown  in 


the  sketch 


We  use  the  notation  given  in  Section  I.  It  is  easily  shown  from  Eq.  (2) 
that  V • u = ^ • ju  =0.  If  we  substitute  Eq.  (4)  into  Eqs.  (1 ) and  (2)  and  use  statistical 
homogeneity  (in  planes)  and  statistical  stationarity  we  find  after  averaging 
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2-1 

J u' v'(z')  dz'  + ^Ml 
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(A-l ) 
(A- 1) 


u'  w'  = v'  w'  =0 
P /2  . . 

f = -ax  - w (z)  (A-2) 

recalling  that  averages  of  terms  in  Eq.  (1)  depend  at  most  on  z.  Here  a is 
the  pressure  gradient  in  the  downstream  direction;  it  overcomes  the  viscous 
drag  of  the  plates  upon  the  flow. 

To  relate  these  results  for  the  parallel  plate  problem  to  the  single 
boundary  layer  problem,  we  replace  the  parameters  a and  h by  the  gradients 
of  the  pressure  and  velocity  at  the  z = 0 plane.  Eqs.  (A-l)  and  (A-2)  be- 
come 

(A- 3) 


a = 


p„  /p 


and 


h = uy  (0)  = -vpu?  «0)/pv  (0) 

d Zj  A 


(A-4) 


where  the  argument  equal  zero  represents  the  value  of  the  function  at  the 
plane  z = 0,  and  the  subscripts  x and  z indicate  the  derivatives  with  respect 
to  those  coordinates.  We  can  substitute  Eqs.  (A-3)  and  (A-4)  into  Eqs.  (A-l) 
and  (A-2)  to  obtain  averages  for  the  related,  simpLe  boundary  layer. 
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Appendix  B 

EVALUATION  OF  SURFACE  FORCES 


We  begin  with  the  stress  on  the  boundary  layer,  exerted  by  the 
wall  on  the  fluid.  The  general  stress  is 


°i3  = P6i3  ' * 

with  the  boundary  condition 


(B- 1 ) 


u = 0,  on  the  wall, 


so 


3u 

dx 


^ V 

— — =0,  on  the  wall. 

ay 


Thus  the  continuity  equation  gives 


3 w _ _ 3u 
" 3x 


— = 0,  on  the  wall 

3y 


From  these  relations  the  force  per  area  exerted  by  the  wall  on  the  fluid  is 


~ H- 


3 u 
dz 


- ' 


a y 
a z 


j + P k 


(B-2) 


where  i,  j and  k are  unit  vectors  in  the  x,  y and  z directions. 

In  the  body  of  the  report  the  image  flow  is  discussed;  it  is  remarked 
that  a force  must  be  added  when  the  wall  is  removed  in  the  imaging  process. 
Thus  a force  term  must  be  added  to  Eq.  (1),  to  represent  the  drag  of  the 
wall;  the  pr  - are  force  is  correctly  given  by  the  pressure  of  the  imaged 
flow.  The  equation  for  the  imaged  flow  problem,  which  replaces  Eq.  (1)  is 
from  Eq.  (B-2) 


u + u 


7u  = - 


_1_ 

P„ 


7p  -i 


V72u 


/ au  ♦ . dv  *1  . 

• ' * Ti  >)  6,zl- 
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with  6(z)  the  Dirac  delta  function.  To  obtain  Eq.  (8)  we  took  the  divergence 
of  the  equation  now  replaced  by  Eq.  (B-3).  Thus  the  wall  stress  gives  an 


*dded  term  to  Eq.  (8)  which  is 

’o  r > / - 1 


u^r')  o2  v'(r^) 


«(z')  dr^ 


L 3z0  Sxo  a«oa*oJ 

Suppose  we  integrate  over  and  compare  estimates  of  the  integrals  of 
(B-4)  andof  Eq.  (8).  Near  the  wall  the  scale  in  the  x and  y directions  is  of 
order  5,  the  thickness  of  the  boundary  layer.  For  (B-4)  the  scale  in  the 
z-  direction  is  the  thickness  of  the  laminar  sublayer  which  is  of  order 

,1/2 


(B-4) 


(B-5) 


Then,  for  example,  the  first  term  in  the  integrand  of  (B-4)  is  of  order, 
after  Zq  integration 


.2  ' 
9 u 


9z' 


TUT 


(B-6) 


6(  V 6/u') 

The  corresponding  quantity  in  Eq.  (8),  after  integration  in  the  Zq  direction,  is 
of  order 

(2 

(B-7) 


,2 

6 T 


or  bigger.  The  ratio  of  (B-6)  to  (B-7)  is  (Re')"1^2 


where 


Re'  = -HlA. 
v 


(B-8) 


If  we  suppose  in  an  application  that  fi  = 15  cm,  u'  = 0.  1 cQ,  with  cQ  the  speed 

2 

of  sound  (the  mean  flow,  tht  a,  near  Mach  1)  and  use  v = 0.  Z cm  /sec,  the 


t- 


-2 

ratio  of  the  terms  is  0.  6 X 10  . The  contribution  of  the  wall  drag  is  leas 

than  1%  of  other  forces  in  the  flow. 

Essentially  the  contribution  of  the  wall  is  a viscous  force,  usually 
(as  here)  of  leas  importance  than  that  of  inertial  forces. 
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Appendix  C 

EVALUATION  OF  A.,  Bt  AND  C 

These  functions  are  needed  for  the  evaluation  of  the  integrals  given  in 
Eq.  (24),  Begin  with  the  functions  A^,  and  refer  to  Eqs.  (26) and  (26a).  A simple 
way  to  find  the  functions  is  the  following.  Take  the  cross  product  of  Eqs.  (26) 
and  (26a)  four  times,  using  the  5_,  and  contract  indices.  Then  the  only  term  in 
Eq.  (26a)  which  doesn’t  vanish  is  the  third  on  the  right  side,  thus  giving  a 
single  equation  for  A^.  Now  continue  by  taking  a cross  product  of  Eqs.  (26)  and 
(26a)  twice  with  5,  summing  on  the  other  two  indices.  This  gives  a single 
equation  involving  A^  and  A^!  substitute  the  value  for  A^  just  found,  and 
determine  A2*  To  determine  Aj,  sum  on  the  four  indices  in  pairs  in  Eqs.  (26) 
and  (26a),  use  the  values  of  A2  and  A^  already  found  and  determine  Aj. 

From  these  manipulations,  one  obtains  the  following  results  tor  these 
func  tions. 

00 

Aj(?)  = S"4  J drr6  Q^r)^  KQ(?,r)  - K2(?,r)  + ^ K4(?,r)j 

0 

00 

A2(5M  T2  / dr  r6  Qf  (r)  [-  | Kq  + | K2  - | (C-l) 

0 L J 

A3(?)  = j dr  .>  Qf  (r)  [jKj-jK2  + ; kJ 
0 

with  (see  below) 

K (§,  r)=  2tt  f X d X_  , (C-2) 

J /a+bx 

a = + r2,  b = -lit 


and 


To  find  the  functions  Ej  and  refer  to  Eqa.  (27) and  (27a),  removing  the 
factor  Using  those  equations  take  the  cross  product  of  both  twice  with 

the  vector  J£,  yielding  a single  equation  for  B.,.  Then  sum  on  the  two  in- 
dices giving  a single  equation  for  Bj  and  B2  from  which  Bj  can  be  found. 
The  result  is 


Bl  = 


B. 


= I"2  f dr  r4  Q1(r)Q2(r)[|-  K2(?,r)  - ^K0(5,r)] 

CD 

Jdrr4Ql(r)Q2(r)[|K0-}  KJ 


(C-3) 


Finally  to  determine  the  function  C refer  to  equation  (29a).  One 
simply  carries  out  the  integral  involved  there  and  obtains  the  result 

CD 

C = J dr  r2  Q2(r)  KQ(5,r)  (C-4) 

0 


The  integrals  in  (C-2)  can  be  found  in  terms  of  elementary  functions: 


K0(5,r) 

K2<S,r) 


5 < r 


4tt 
15  r 


5+2 
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§ < r 


K4(S,r) 


4tt 

tut 


63  + 36 


5 < r 


L r r j 

and  are  eymmetric  in  § and  r , as  must  be  so  from  their  definition,  (C-2). 
The  functions  have  discontinuous  derivatives,  which  lead  (under  two  differ 
entiations)  to  Dirac  delta  functions  in  the  work. 
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Appendix  D 


INTEGRATION  BY  PARTS,  AT  A SINGULAR  POINT 


It  will  be  convenient  to  integrate  Eq.  (8)  by  parts  (i.  e.  , use  Gauss' 
theorem).  In  the  process,  the  singular  point  at  r*  - r^  = 0 must  be  treated 
carefully;  we  present  the  reasoning  in  detail.  First  note  that  as  it  stands, 
the  integrand  is  integrable  at  that  point.  Thus  we  may  set  up  a small  sphere 
of  radius  c about  the  point  and  delete  the  integral  over  sphere- e from  the 
rest  of  the  volume  of  integration  (all  space)  without  changing  the  result.  Let 
the  remaining  volume  be  V-,  whereas  all  space  we  let  be  V.  Thus  we  have 
for  Eq.  (8) 


p'(i'»*5r  JJ/lio  -xT1  8 ;■  n»e(iold^ 

t t a.  no 
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Use  Gauss'  theorem  on  the  first;  we  have  two  surface  integrals,  one  at 

infinity  (i.  e.  , large  distance)  which  vanishes  because  all  fluctuations  vanish 

there,  and  the  other  on  the  e-sphere;  that  one  vanishes  as  well  because  the 

2 

element  of  surface  area  is  of  order  e , giving  an  integrand  of  order  e.  Thus 
the  first  term  on  the  right  side  of  (D-2)  vanishes.  Consider  the  second: 
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Use  Gauss'  theorem  on  the  first  term  on  the  right  side;  it  is  equal  to  (recalling 

that  surface  integrals  at  infinity  vanish) 

Po 


£ Jl  «ws  h„b<£o>  TT- 1— o - il ''  ds<4>  <D-4> 

€ ^0* 

Let  £ = £.0  " — r an<*  reca^  that  f|  is  directed  into  the  surface  of  the  sphere 

and  find  for  (D-4)  for  c very  small: 
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with  dO  the  element  of  solid  angle.  This  is  easily  seen  to  give 
Po  4tt  / 

— 4n*  — **aQ,(£  )•  Combining  these  results  we  find 
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Consider  now  the  value  of  the  integral  in  Eq.  (D-6),  if  carried  out 
within  £.  We  see 
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using  the  vector  properties  of  the  integral.  Evidently  if  we  (incorrectly) 
drop  the  surface  integral  which  gave  the  added  term  to  Eq.  (D-6)and 
(incorrectly)  extend  the  volume  integral  in  Eq.  (D-6)  throughout  e,  we  have 
cancelling  errors  and  the  correct  result 
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(D-8) 


Some  care  must  be  exercised  to  be  sure  that  the  delta -function  contribution 
at  = j/  is  not  lost. 


-61- 


Appendix  E 

TECHNIQUES  USED  IN  EVALUATION 
OF  THE  PRESSURE  CORRELATION 


In  order  to  numerically  evaluate  the  pressure  correlations 
Eqs.  (38)  and  (44)  the  angular  integration  is  first  performed  analytically 
in  terms  of  Hypergeometric  functions.  From  Eq.  (38) 


iiiLifi 

|r-5|3 


a"'5  . 75|(^-r2)2b'2  A(x)  + 2(52-r2)b"'  B(x)  + C(x)|  -52b_1B(x) 


(E-l) 


where 


x = ^ , a = (p  + T])2  + (z  - 53)2,  b = (p-n)2  + (Z  - ?3)2 


and 


C(x)  = 2 

B(x)  = 2 


(1-y2)"* 5(l-xy2)‘* 5dy  = 2 K(x) 
(1-y2)  * 5(l-xy2)"*  5dy  = 2E(x) 


A(x)  = 2 [(2 -x)B(x)  5(l-x)C(x)]/3  (E-2) 

Note  that  ?2  = p2  + ?2  and  r2  = r\Z  + z2. 

O 

Simple  polynomial  expressions  accurate  to  two  parts  in  10°  exist  for  the 
complete  elliptic  integrals10  K and  E.  Note  that  for  x = 0 (for  example) 
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calculation  of  the  variance  or  the  correlation  for  z separation  only), 

2 

A = B = C = tt.  For  r equal  to  0 (i.  e.  the  variance)  Eq.  (E-l)  becomes 
simply  2n/5.  This  implies  that  the  fourth -order  variance  is  given 
simply  by 


PP 


o 

2TT 


0© 

-[u/(z/)]4D2(o)  + f d5_[u'(  |z  ' + 5,  | ]4M(xx)m"1 


(E-3) 


where 


M(xx) 


dSD^?) 


(E-4) 


and 

xx  = §3M_1 

M = (T"1  + |z'  + ?3  I “ 1)“ 1 

The  function  M(xx)  is  a universal  curve  independent  of  the  velocity  field  and 
should  not  be  confused  with  the  scale  factor  M.  Continue  with  the  calcu- 
lation of  the  correlation  C^5)  W substitution  of  (E-l)  into  Eq.  (38).  Direct 
numerical  integration  at  this  point  leads  to  significant  error  and  large 
computation  times  due  to  the  apparent  singularity  near  £ = jr . Instead 
change  integration  variables  letting 


Expression  (E-l)  becomes 


[(n2-»2)(A2-a>2)-4rizAu)]  B t 2(,n2A2-2T|zAu)-fu)  2z2  [B(x) 5C(x)] 

b4  ab 

+ [3A(x)  + B(x)]  + [.  75A(x)  + . 5 B(x)  + . 75  C(x)] 

b 


Substitute  (E-5)  into  Eq.  (38)  and  integrate  the  first  term  in  (E-5) 
by  parts  w.  r.  t.  A.  This  gives 
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V (z 


*1.  z)  = 
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00 
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dujfu^  |z' | )]' 


M 


00 

/ 

-T1 


dA  P(r),  z,  z'  A,  id) 


where 


Z»  Zq»  A,  Ul)  - 


[ [(T12-z2)A-2tlzU)]  f -.5 


-1.  -1, 


J^a  BDj  + 2a  * pTpc-  (l-2p(p+T))a 


X (B -C)Dj  + p ? I zq  I a-*  BDj  - pfp+TDa"1,  5BDjJ 

+ [l!3Al±SJ  x 2(t12A2-2t)zAuj+U)2z2)  (B-.5C1 

L b ab 

+ (.  75A  + . 5B  +.  75C)  J(a'*  5pDj) 


(E-5) 


(E-6) 


(E-7) 


a = (211  + A)2  + uu2 
b = A2  + ii>2 
§ = (P2  + §2)‘5 


% = * " 


U) 


x = 4TKri  + A)a 


-1 


Xq  - X*  + UJ 

The  argument  of  the  functions  A,  B,  and  C is  x while  the  argument 

of  Dj  and  the  derivative  of  "Dj"  is  ?M  To  complete  the  numerical 

evaluation  of  Eq.  (E-6)  the  integral  over  A was  first  broken  with  two 

terms,  one  from  -r\  to  T|  and  the  other  from  r)  to  ®.  Then  the  integral 

over  A from  -r|  to  f|  was  rewritten  in  terms  of  a single  integral  from  0 to 

Ti.  Finally  the  integral  over  U)  from  -®  to  ® was  rewritten  as  a single 

integral  from  0 to  ®.  Note  that  [u#  (|  z^|  )]^M~2  was  brought  inside  the 

A integral  in  the  final  expression.  The  final  results  do  not  give  numerical 

difficulty  at  A = uu  = 0 and  typical  evaluation  times  for  T),  z not  too  large  are 

of  the  order  of  10  seconds  on  a CDC  7600.  Note  that  the  integral  over 

A is  a universal  function  of  (T|/M,  z/M,  id/M).  If  one  knew  this  function, 

(2) 

it  would  be  very  easy  to  evaluate  ^pp  for  any  velocity  field.  Unfortunately 
this  function  is  too  complicated  to  compute  and  store  for  all  parameter 
values.  Instead  one  may  consider  the  special  cases  when  either  r|  or 
z are  equal  to  zero.  This  results  in  a two  parameter  function  which  can  be 
tabulated  and  stored.  If  there  are  many  different  boundary  layer  cases  to 
be  considered,  there  may  be  a significant  advantage  of  tabulating  these 
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functions.  However,  even  on  a CDC  7600,  the  calculation  time  required 

2 1 

to  determine  these  functions  for  C and  C terms  will  take  approximately 

PP  PP 

one  hour.  Proceeding  now  with  the  case  when  T)  = 0 the  fourth-order 
correlation  becomes 
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(E-8a) 


where  from  Eq.  (E-7) 


and  a change  of  variables  was  made  within  Eq.  (E-8b)  so  that 
b = A2  + KZ 

% = ^A2  + (B-A)2 
A = U)/M 
B = z/M 


(E-8b) 


In  obtaining  Eq.  (E-8)  the  partial  integration  proceeded  in  a slightly 


different  way  than  discussed  after  Eq.  (E-7)  since  a = b.  Depending  on 
the  accuracy  required,  this  function  Mean  be  stored  in  tabular  form 


srrTwr  R'win  "n ,T"  yj  #*  * ^nwrii 


or  the  complete  double  integral  Eq,  (E-8a)  can  be  calculated.  If  the 
first  method  is  used  an  interpolation  is  performed  for  intermediate 
values  of  A and  B . In  a similar  manner,  for  z = 0 
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where 


00 

M(A,  0,  C)  = M-1  / dAp  (n,  0,  U'|,  A,  id)  = 


-T) 


/( 


^Ab'V*5 


B(x)D1(5)  + 2a"1x"1pC[l-2p(P+C)a_1][B(x)-C(x)]D1(5) 


+ p2?_1B(x)Dj(§)  - p(p  + CJa^BfxJD^S) 


C Ab-1[3A(x)  + B(x)]  + 2C2A2b'1a"1  [B(x)  - . 5C(x)]  + 


a_1[.  75 A(x)  + . 5B(x)  + . 75C(x)] 
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and 
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5 

a 
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u)/M 
11/M 
A2  + A2 
A + C 


^TTs2 


(p  + cr  + a 

4pC  a" 1 


2 


Again  the  integration  over  A should  be  broken  into  two  parts,  one  from  0 
to  Cand  the  other  from  C to  “ for  use  of  numerical  integration. 

The  calculation  of  proceeds  in  a similar  way  though  it  is 
much  easier.  First  the  integration  is  performed  over  angle.  The 
required  scalar  integrals  are  [the  underbars  do  not  represent  vectors 
on  A (x)  and  B (x)] 


a'1/2C(x)  = f jr-il"1 
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TT 

a"l^2A(x)  = a"1^2[C(x)-4B(x)]  = j cos20|r-?|  ~l 


d0 
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Now  B(x)  is  not  easily  expressed  in  terms  of  the  elliptic  integrals; 
therefore,  B(x)  is  calculated  as  a function  of  x for  04xM  and  stored 
for  later  use.  A polynomial  fit  could  also  be  made.  Substituting  Eqs. 
(E-10)  into  Eqns,  (42)-(44)  yields 


^ 2 00  00 

cin(ri'0'z,z,)  = ~"n“  pdp  N(p,U)  ,r),  z,  z' ) 
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and  it  is  recalled  that 
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Replace  A (x)  by  B(x)  in  C^}  (0»  Tl»  z,  z').  Direct  numerical  integration 
of  Eq.  (E-l  1)  may  follow  though  there  may  be  a time  advantage  by  repla- 
cing the  integral  over  id  from  -»  to  00  by  one  from  0 to  “>.  When  either  z, 
or  r]  is  zero,  then  can  be  determined  in  terms  of  a two  parameter 
scalar  function  N which  is  equal  to 


N(AT,  B,  C)=  M"j  pdp  N(p,  a),  n,  z,  | z'\  ) = 

0 0 

00 

M-Zpdpfd-^  j C7DB(?)  - Q1(5>]  p2r2A(x)  - 

DB(§)C<x)|  + [f(?)  - Q2<?)  - DB(  ?)]  P2?'ZA(x) 

- t2f(?)  - ?2DB(?)]  C(x)/5^a"*  5 

where 

A = u'/M 

B - z/M 

C = r|/M 

a = (P  + C )2  + A 2 

5 = '/"p"2  + (B  - A)2 

x = 4 pC  a"1 


(E-13) 


and  A (x)  and  B (x)  are  interchangeable  depending  on  whether  the  x or  y 
correlation  function  is  desired.  The  variance  is  easy  to  determine  since 
then  N reduces  to  a one  parameter  function. 


Appendix  F 


GRAPHICAL  REPRESENTATIONS  OF  VARIOUS  FUNCTIONS 
USED  TO  CALCULATE  THE  PRESSURE  CORRELATION 

Figure  F-l  is  the  plot  of  the  longitudinal  velocity  coefficient 

f(r)  vs  r on  a log-log  scale.  Figure  F-2  is  the  plot  of  g(r)  vs  r on  a 

linear  scale.  A plot  of  Q^(r)  vs  r is  shown  on  a log-log  scale  in  figure 

F-3.  ^{r)  and  -D^(r)  vs  r are  shown  in  figs.  F-4  and  F-5  using  log- 

-7— 

log  scaling.  The  quantity  5 E^S)  vs  5 is  shown  in  figure  F-6  on  a 
log- log  scale.  The  function  B (x)  vs  x is  shown  in  figure  F-7  on  a linear 
scale.  Note  that  tabular  data  is  included  in  figures  F-l  through  F-7. 

As  mentioned  in  Appendix  E,  it  is  possible  to  obtain  two  para- 
meter functions  Mand  N,  when  either  T]  or  z is  equal  to  zero,  from 
which  the  correlation  functions  could  be  calculated.  For  a large  number 
of  correlation  calculations  for  differing  flow  field  conditions  this  technique 
would  have  a significant  computer  time  advantage  over  performing  a 
double  integral  for  each  desired  point.  For  T]  - 0 the  function  Kl  (A,  B,0) 
is  calculated  and  shown  in  perspective  in  figure  F-8.  The  parameters 
A = u)/M(x  axis)  and  B = z/M  (y  axis)  with  M = (T  * + | z/  + <u| '*)"’*  with 
ranges  of  -100  £ A £i  and  Os  B £ 100  were  used  to  calculate  M.  The 
parameter  interval  was  varied  logarithmic ly.  Figure  F-8  views  the  function 
from  a location  of  x = 10,  y = 100,  z = 100.  The  plot  shown  is  for  the  range 
of  values  -5  s A si  and  0 s B s 5.  The  scaling  is  reduced  by  a factor  of  5 
in  the  y direction  and  a factor  of  10  in  the  z direction.  Note  that  there  are 
two  sheets  to  the  function  - one  for  positive  A and  the  other  for  negatived. 
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Fig.  F-5.  D'  j (r)  vs  r 
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The  values  for  X = 0 have  been  removed  to  show  a clean  figure.  If  KT  is 

|2U\ 

used  to  calculate  C'  ' the  two  sheets  are  calculated  separately  and  then 

PP 

added  together  during  the  computation.  Figure  F-9  is  the  positive  sheet 
of  W (X",  B,  0)  alone  from  the  vantage  point  of  x = 10,  y = 10,  z - 10. 

The  z scale  has  been  reduced  by  a factor  of  10  and  W is  only  plotted  for 

0 £ S £ 1 and  0 s B si.  Figure  F-10  contains  both  positive  and  negative 
bheets  of  M (X,  B,  0)  plotted  separately  from  the  vantage  point  x = -10, 
y = 10,  z = 100  and  a z scaling  reduced  by  a factor  of  10.  From  figures 
F-8,  F-9,  and  F-10  it  is  seen  that  the  S positive  sheet  is  maximum  at  the 
origin  and  falls  rapidly  for  increasing  values  of  S’  or  B.  The  S negative 
sheet  on  the  other  hand  increases  with  increasing  B , reaching  a maximum 
value  near  B = . 1 then  falls  rapidly  to  zero.  The  other  available  fourth- 
order  contribution  is  for  z = 0.  Then^T  (X  , 0,  S ) is  required  where  C = 
ti/M.  The  function  was  calculated  over  the  range  0 sXs  10  and  0 sSs  10. 
Figure  F-ll  is  a perspective  plot  of  this  function  observed  from  a position  of 
x = 100,  y = 100,  z = 100.  Figure  F-12  is  the  same  function  observed  from 

x = 100,  y = 100,  but  z = 10.  It  is  seen  that  the  function  is  highly  peaked  at 

X = <5=  0 falling  rapidly  as  both  X and  C increase.  The  function  falls  a bit 
faster  in  the  y(C)  direction  becoming  negative  then  approaching  z = 0 from 
the  negative  values.  In  the  special  case  that  p or  z equals  zero  N (Eq.  E-13) 
can  be  determined.  Figure  F-13  is  a perspective  plot  of  V (X  , B , 0)  (t|  = 0) 
viewed  from  from  a location  of  x ® -20,  y = 10. , z = .3  with  the  z scaling 
increased  by  a factor  of  100.  The  range  of  values  are  for  -10  <x(X)  £ 1 and 
0 * y(B)  <■  10.  The  same  function  is  viewed  from  a location  x = -20,  y = -20, 
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Fig.  F-13.  Perspective  of  5T(A,  5",  0)  Viewed  from  a 
Location  of  (-20,  10,  .3).  The  z scaling 
has  been  increased  by  a factor  of  100. 
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z a 0 and  x = 0,  y = 20,  z = . 3 in  figures  F-14  and  F-15  respectively. 
Perspective  views  of  N (R  , 0,  C ) are  given  for  the  x correlation  function  in 

figures  F-16  and  F-17  where  the  viewing  locations  are  x = 20,  y = -20,  z = 

. 3 and  x = -20,  y = -20,  z = . 3 respectively.  The  z scaling  has  been 
increased  by  a factor  of  50.  The  perspective  views  of  N (A,  0,  75 ) for  the 
y correlation  function  are  giver,  in  figures  F-18  and  F-19.  In  figures  F-16  - 
F-19  the  figures  are  plotted  for  0 £ x(5)  £ 10  and  0 £ y(C)  £ 10. 
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Fig.  F-17.  Perspective  Plot  of  N(A,  0,  C) 
(x  correlation)  Viewed  from  a 
location  of  x = -20,  y - 20, 
z = 0,  3.  z scaling  increased 
by  a factor  of  50. 


Fig.  F -18.  Perspective  View  of  ?T(A,  0,  (T)  for  the 
y Correlation  Viewed  from  the  Location 
x = 10,  y = -20,  z = 0.  3.  The  z scaling 
is  increased  by  a factor  of  50. 


Fig.  F-19.  Perspective  View  of  IT(X,  0,  Cl  for 
the  y Correlation  Viewed  from  a 
Location  of  x = -20,  y = -20,  z-  0.  3. 
The  z scaling  is  increased  by  a fac- 
tor of  50. 


THE  IVAN  A.  GETTING  LABORATORIES 

Tha  Uboritory  Operation*  ol  Tha  Atrosptct  Corporation  ia  conducting 
experimental  and  theoretical  Investigation*  naceaaary  lor  tha  (valuation  and 
application  of  aeiantiflc  advaneaa  to  naw  military  concept*  and  systems.  Var- 
aatility  and  flexibility  have  bean  developed  to  a high  degree  by  the  laboratory 
personnel  ia  dealing  with  tha  many  problem*  ancountarad  in  tha  nation'*  rapidly 
developing  *pec*  and  miialla  ayttam*.  Exparti**  in  the  latest  aeiantiflc  devel- 
opment* ia  vital  to  tha  accomplishment  of  teak*  related  to  tha**  problem*.  Tha 
laboratories  that  contribute  to  this  rsaaarch  ai a: 

Aaroohvlc*  Laboratory)  Launch  and  reentry  aerodynamic*,  heat  t ran* - 
t r,  reentry  physic*.  chemical  kinetic*,  structural  mechanic*,  flight  dynamic*, 
a oioepharic  pollution,  and  high-power  gas  laser*. 

Chemistry  and  Physic*  Laboratoryi  Atmospheric  reaction*  and  atmos- 
pheric  optic*,  chemical  reaction*  in  polluted  atmosphere*,  chemical  reaction* 
of  excited  specie*  in  rocket  plume*,  chemical  tharmodynamica,  plasma  and 
laser-induced  reaction*,  lasar  chomistry,  propulsion  chemistry,  space  vacuum 
and  radiation  affacta  on  matarial*,  lubrication  and  surface  phenomena,  photo- 
sensitive  material*  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  phyalca  and  chomistry  to  problem*  of  law  anforcamant  and  biomedicine. 

Electronic*  Research  Laboratoryi  Electromagnetic  theory,  device*,  and 
propagation  phenomena,  including  plasma  electromagnetic*;  quantum  electronic*, 
lasar*.  and  electro-optic*;  comm-wication  aclenc**,  applied  electroale*,  semi- 
conducting, superconducting,  and  crystal  device  phycica,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Science*  Laboratoryi  Development  of  now  materials;  metal 
matrix  composite*  and  saw  forma  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  la  reentry;  spacecraft  material*  and  electronic  component*  in 
nuclear  weapon*  environment;  application  of  fracture  mechanics  to  *tree*  cor- 
rosion and  fatigue -‘nduced  fractures  in  structural  metals. 

Rose*  Rcloacs*  Laboratory:  Atmospheric  and  ionospheric  physic*,  radia- 
tion from  the  atmosphere,  deadly  and  composition  of  die  atmosphere,  aurora* 
and  airglow;  magaetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnet  ; sphere;  solar  physics,  studies  of  solar  magnetic 
A elds:  apace  astronomy,  x-ray  astronomy;  th*  affect*  of  nuclear  axploeicn*, 
magnetic  (term*,  and  solar  activity  on  the  sarth's  atmosphere,  ionosphere,  and 
magnetosphere:  the  affect*  of  optical,  electromagnetic,  and  particulate  radio- 
Won*  In  spec*  on  space  system*. 
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